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Abstract 

The interface problem for the linear Korteweg-de Vries (KdV) equation in one-dimensional 
piecewise homogeneous domains is examined by constructing an explicit solution in each domain. 
The location of the interface is known and a number of compatibility conditions at the boundary 
are imposed. We provide an explicit characterization of sufficient interface conditions for the 
construction of a solution using Fokas’s Unified Transform Method. The problem and the 
method considered here extend that of earlier papers to problems with more than two spatial 
derivatives. 


1 Introduction 

Interface problems for partial differential equations (PDEs) are initial boundary value problems for 
which the solution of an equation in one domain prescribes boundary conditions for the equations 
in adjacent domains. In applications, interface conditions are often obtained from conservation 
laws m- Few interface problems allow for an explicit closed-form solution using classical solution 
methods. Using the Fokas or Unified Transform Method HI El IS], such solutions may be constructed 
for both dissipative and dispersive linear interface problems as shown in 0 0 H2J HU H3 HU US]. 

All previous papers addressing interface problems using the Fokas Method for interface problems 
have dealt exclusively with problems that are of second order in the spatial variable. This paper is 
the first investigation into higher-order problems. The process presented in this paper makes clear 
how to resolve new issues that arise when moving to a higher-order problem. 

The nondimensionalized Korteweg-de Vries (KdV) equation 

Qt T QqQx T Qxxx — 0, 

is one of the most studied nonlinear PDEs mmmm- It arises in the study of long waves in 
shallow water, ion-acoustic waves in plasmas, and in general, describes the slow evolution of long 
waves in dispersive media [2]. In this manuscript we study the linearized KdV equation (LKdV) in 
a composite medium, 

Qt = cr(x)q xxx , - oo < x < oo, (1) 

where cr(x), a real-valued function, is piecewise constant. This equation describes the behavior 
of solutions of the KdV equation in the small-amplitude limit and understanding its dynamics is 
fundamental in understanding the dynamics of the more complicated nonlinear problem. 


1 


In what follows an explicit solution method is given resulting in closed-form expressions. We 
provide criteria which, under the additional assumption of existence of a solution, are sufficient 
for a solution representation to be obtained via the Fokas Method. Although we do not prove 
uniqueness of the solution, we note some examples of interface conditions that do yield uniqueness. 
The numerical evaluation of the solution is not considered but should be possible via the methods 
presented in mmmm- As we do not have a physical application on hand, this paper addresses 
the mathematical question of the number and type of interface conditions required to ensure that ([!]) 
is well posed. 

2 Background 

Determining the number of boundary conditions necessary for a well-posed problem is a nontrivial 
issue, especially for boundary value problems (BVPs) with higher than second-order derivatives. 
Consider LKdV posed on the half line 


qt = cr 3 Qxxx, x>0, t > 0, 


( 2 ) 


where the form of the coefficient a 3 is chosen for convenience. If a < 0 then one boundary condition 
is needed, whereas if a > 0, two boundary conditions must be prescribed in order for the problem 
to be well posed m- This difference in seemingly very similar BVPs is understood at an intuitive 
level by considering the phase velocity c(k) = — iu(k)/k where a ){k) = icr 3 k 3 [TT] , Thus, the phase 
velocity is c(k ) = a 3 k 2 . If a < 0 the phase velocity is negative and information travels toward the 
boundary as in Figure [Taj If a > 0, the phase velocity is positive and information travels away 
from the boundary as in Figure lb Therefore, it seems reasonable that one must prescribe more 


boundary information. Note that if we were solving ([2]) for x < 0 these results would be switched. 
This will become relevant in what follows for the interface problem on the whole line. 


t t 




Figure 1: (a) When a < 0 information from the initial condition propagates toward the boundary 
x = 0 and one boundary condition needs to be prescribed, (b) When a > 0 information from 
the initial condition propagates away from the boundary and two boundary conditions need to be 
prescribed. 
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Remark. The above argument hints at the Method of Characteristics. This analogy is not 
justified, as the problem at hand is dispersive and energy spreads along phase or group velocity 
rays, rather than travels along them. As such, the above argument is entirely heuristic but, with 
hindsight, it provides some intuition. 


One of the strengths of the Fokas Method for solving linear PDEs is the straightforward way 
it enables determination of the number and type of boundary conditions required for a well-posed 
problem 13 Eli. Previous papers by us and others & Ema m m na m have shown that 
the Fokas Method is useful for finding explicit general solutions to interface problems. In the cases 
currently in the literature, only second-order problems are considered and the number of conditions 
required at each interface is clearly two. The example of LKdV on the half-line suggests that the 
number of interface conditions needed in the case of LKdV with an interface depends on the sign 
of a. This is the case indeed. In Propositions 0® we describe exactly the number and type of 
conditions necessary. 


3 Notation and Set-Up 


We investigate 0 where cr(x ) is the piecewise constant real-valued function 


a(x) 


af, x < 0, 
erf, X > 0, 


(3) 


with the initial condition q(x,0) = qo(x), and appropriate conditions at the interface x = 0. The 
choice of the power 3 in the definition of a(x ) is purely for convenience. We assume throughout 
this work that the solution decays rapidly to zero as |x| —> oo. If nonzero conditions at |x| = oo 
are desired this can be treated easily in a manner similar to that for the heat equation in [BJ and 
for the linear Schrodinger equation in m- We pose ([!]) as the following interface problem: 


Qt ^lQxxxi a? < 0, 

„(2) 3 (2) > n 

Ht ( ^2 ( ixxxf X U, 

subject to the initial conditions 

Q (1 \x, 0) = gj 1} (x), 
q {2) (x, 0) = q ( 0 2 \x), 


0 < t <T, (4a) 

0 < t < T, (4b) 

x < 0, (5a) 

x > 0, (5b) 


with gW(-, t) G 5( —oo,0) and G 5(0, oo) where S(X ) is the Schwartz space of restrictions 

to X of rapidly decaying functions. Likewise, we assume rapid decay of the initial conditions, 
(/q 1 )(•) G S(— oo,0) and q^ 2 ' (■) G 5(0, oo). Note that we do not require (0) = q { 2 ' 1 (0), rather we 
assume that the initial data are compatible with the interface conditions, which are specified below. 

Some number of interface conditions at x = 0 needs to be prescribed. The number and type of 
such conditions are given in Propositions [l}|3| We make a distinction in this manuscript between 
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“boundary problems” and “interface problems.” Boundary problems are those in which the condi¬ 
tions given at the interface x = 0 allow one to solve either (4a) or (4b) as a half-line BVP without 
knowing the solution on the other domain. For example, if one can solve a BVP for q^\x, t) then 
one can use that solution to provide any necessary conditions at x = 0 to solve the second BVP 
for q^ 2 \x,t). Conditions for a well-posed BVP are given in 0 (22]. Since these cases have been 
examined, we restrict to those interface conditions which do not decouple such that either (4a) 
or (4b) can be solved as a BVP. 

It is of note that by making use of the PDE, interface conditions can always be written as 


linear combination of 


d n 
dx r ‘ 


Z (1) (M) 


x=0 


d n 

and —— q^ 2 \x,t) 
dx nH v ’ ’ 


for n = 0,1, 2 


x=0 


= f(t), 


for all t. For example, one might require q^J x (0, t) = q^JxiOX) as an interface condition. This can 
be imposed by applying the equation and integrating in t to give 


v 

*1 


(o,t) 


- 4« ,2) (o.t)= p 9 P(o) - 44 2, (o), 


<To 






( 6 ) 


for all t, which is clearly of the form we require with f(t ) = (0) — \q^P (0). Using a similar 

process for any conditions on derivatives greater than second order as well as elementary linear 
algebra one can always express the interface conditions in the reduced forms given in Propositions[l]- 
[3] possibly after letting x —»• —x. 


Remark. If an interface condition specifies a linear combination of d^q^(0,t) (n=0,l,2) only 
or <9^ {2) (0 ,t) (n=0,l,2) only, then we say it is a boundary condition. Note that the interface 
conditions 

q( 1 \0,t)=0, and q^{0,t) — q^ 2 \0,t) = 0, 
are equivalent to the interface conditions 

g^(0, t) = 0, and q^ 2 \0,t)=0, 

so it is only meaningful to discuss the maximum number of boundary conditions for any equivalent 
expression of a given system of interface conditions. Henceforth any mention of a number of 
boundary conditions should be interpreted as such a maximum number of boundary conditions. 

A problem with one boundary condition may or may not decouple into a pair of BVP. Even if 
such a decoupling is possible, it may or may not be possible to solve the BVPs sequentially. For 
example, the problem with (j\. a 2 > 0, boundary condition q^\0,t) = 0, and interface conditions 
q x \o, t) = q x 2 \o,t) and f/ix (0, t) = q x 2 J (O.t) decouples into a solvable BVP for q ^ and a subse¬ 
quent solvable BVP for q^ 2 ' 1 . However, the problem with <7 i,< 72 < 0, and the same boundary and 
interface conditions does not decouple. 
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Figure 2: Regions for the application of Green’s Formula in the case of two semi-infinite domains. 


4 Application of the Fokas Method 

We follow the standard steps in the application of the Fokas Method. Assuming existence of a 
solution, we begin with the so-called “local relations”: 

{e - ikx+UJj t q U) )t = (e -to+o,, V 3 {q U) + ikq U) _ k 2 q{ i) ))xj 

with Uj = Uj(k ) = ia'jk 3 for j = 1, 2. Applying Green’s Theorem [T] and integrating over the strips 
(—oo,0) x (0,f) and (0,oo) x (0, t) respectively (see Figure [2]) we have the global relations 


0 -ikx+^t q { i) ( X) dx = f e -ikx g ( i) dx 


' —OO 


+ / e"'V? (gWfO.s)+!ti 1) (0,s) - 1», 


r 

Jo 


roc 

3 -ik X +U 2 t q ( 2) dx = / e -^ g J 2 ) ( x ) dx 

Jo 

t 


- / e" 2 Vf (0, s) + i/c4 2) (0, s) - k 2 q {2) (0, s)) ds. 


Let C + = {z € C : 1111 ( 0 ) > 0}. Similarly, let C = {z £ C : Im(z) < 0}. Define k = kn + ikj 
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where kn, kj E M. We define the following: 


II 

/ e~ ikx q W {x,t)dx, 

J — OO 

kj > 0, 

0 < i < T, 

?o 1} (*0 = J 

/*0 

1 e~ ikx q ( o\x) dx, 

kj > 0, 


q {2 \k,t) = 

/‘OO 

/ e- ifc 'V 2) (*,i)dx, 

Jo 

kj < 0, 

0 < i < T, 

Qo\k) = 

roc 

/ (x)dx, 

0 

kj < 0, 


II 

'•+ 0 ' 

e 


n = 0,1, 2, 

0 < i < T, 

II 

-+o 

e 


n = 0,1, 2, 

0 < t < T, 


The global relations become 

e Ult q^(k,t) = q^\k) + a\ (g 2 (uji,t) +ikgi(ui,t) - k 2 g 0 (uji,t )) , ki > 0, (7a) 

e U2t q( 2 \k,t) = q<?\k) - a 3 (h 2 (u 2 ,t) + ikhi{u 2 ,t) - k 2 h 0 (uj 2 ,t )) , ki < 0. (7b) 

We wish to transform the global relations so that g n {-,t ) and h n (-,t ) for n = 0,1, 2 depend on a 

common argument, ik 3 as in [31 E], Noting ?'/c 3 is invariant under the transformations k —> ak and 
k —> a 2 k where a = e 2mt / 3 anc [ evaluating at t = T we have the following six global relations: 


02 


e lkiT q {1) + { a i 92 (ik 3 , T) + ia J kafgi(ik 3 , T) - a 2j k 2 aig 0 (ik 3 , T)) , 

ci Im^'fc) > 0, 

gife 3 T^(2) = ~(2) _ ( a %h 2 (ik 3 ,T) +ia J ka 2 h 1 (ik 3 ,T) - a 2j k 2 a 2 h 0 (ik 3 ,T)) , 

c 2 Iiri^/c) < 0, 


for j = 0,1,2. 

Inverting the Fourier transform in Q we have the solution formulas 
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*qo (k) dk 


e ikx—u>it~(l) i 

—oo 

c 3 r 00 

+ 2tt e lkx ~ uit (g 2 (ui,t)+ikgi(ui,t) - k 2 go(ui,t)) dk, 


(9a) 


1 /*°° 

g (2) (x, t) = — e ikx ~ U2t %’ (k) dk 
J_ oo 


_ ^2 / glfca; w 2 t (h 2 (u 2 ,t) + ikhi(u 2 ,t) - k 2 ho(uj 2 ,t)) dk, 


— OO 

3 /*oo 


(9b) 


2vr 


6 




for t E (0, T) and x E (— 00 , 0) and x E (0, 00 ) respectively. Next, we transform the second integral 
in each of the previous equations so that g n (-,t) and h n (-,t) depend on ik 3 for n = 0,1, 2. 




0 ikx—u)\t ~( 1 ) 




i—x—ik 3 t 

e 47 1 


<7o ' ( fe ) 

(a\g 2 {ik 3 ,t) + ikaigi(ik 3 , t) — k 2 go(ik 3 ,t )) dfc, 


(10a) 




Jkx-u 2 t ~( 2 )( fc ) 


1 

27T 


% 


i — x—ik 3 t 

e CT 2 


{p2h2(ik 3 ,t) + ik<T2hi(ik 3 , t) — k 2 ho(ik 3 ,t )) d/c, 


(10b) 


Let h = {fc£C: Re(ifc 3 ) < 0} = L>tf) U D® U L>( 5 ) as in Figure [3j The parenthetical numbers 
in the superscript of D represent an enumeration of the sectors of the complex plane, in contrast to 
the parenthetical numbers in the superscript of q (and T, below), which represent the two half-line 
domains (— 00 , 0) and (0, 00 ). Let 


D r = {k E C : k € D and \k\ > R} = D& U D^ U D 


)(3) | , n (5) 


'R ’ 


where R > 0 is a positive constant as shown in Figure 0 Let rtf) be the contour d{k E Dr : 
(—1 ) J Ojki > 0}, oriented so that hi 1 ' and lie to the left, and lies to the right of any rtf) 


to which they are adjacent. Note that whether rtf) is the boundary of hi 1 ' Uhl 3 ' or hn' depends 


not only on j but also upon the sign of cr ? . The integrand of the second integral in (10a) is analytic 
and decays as k —> 00 from within the set bounded between M and T^ 1 ), and the integrand of the 
second integral in (10b) is analytic and decays as k —> 00 from within the set bounded between M 


and rtf). Hence, by Jordan’s Lemma and Cauchy’s Theorem, the contours of integration can be 
deformed from M to rtf). 


1 r°° 

q{1) ( x , t ) = -j e ikx -^\k)dk 


1 f 

+ —- / e ff i 
27T J r (i) 


i—x—ikPt. 


{a\g 2 (ik 3 , t) + ika\gi(ik 3 , t) — k 2 go(ik 3 , t)) d k, 


I / ou 

7 ( 2 ) ( Sj t) = — e ikx -^ tq g ] (k) d k 


2vr 


(11a) 


(lib) 


-— [ e* CT 2 X lk * (a 2 h2(ik 3 ,t) + ika2hi(ik 3 ,t) — k 2 ho(ik 3 ,t)) dk. 

27r J T ( 2 ) 

We replace t by T in the arguments of gj and hj by noting that this is equivalent to replacing the 
integral f* e lk3 -j^q^ (0, s) ds with / Q T e lk? ‘ (0, s) d s-ff e* fc3 J^tf) (0, s) ds. Using analyticity 

properties of the integrand and Jordan’s Lemma, the contribution from the second integral is zero 
and thus, 




1 /*°° 

— / e ikx - Wlt q { o\k)dk 
2vr J_ 00 

+ [ e^ x ~ ik3t (aj 9 2(ik 3 ,T) + ika m (ik 3 ,T) - k 2 g 0 (ik 3 ,T )) d k, 

27r J r (i) 


(12a) 
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Figure 3: The evenly distributed regions D^ where Re(ik 3 ) < 




Figure 4: The regions D^\ D^ \ where Re(*A: 3 ) < 0 and \k\ > R. 










(12b) 


7 ( 2 ) 


i r 

{x ' t)= ^L 


e ikx ~ U2t q^)( k ) d k 


i—x—ik 3 t 

e CT 2 


(a^hhik 3 , T) + ika 2 h\(ik 3 , T) — /c 2 /io(ifc 3 , T)) dk. 

2vr J r ( 2 ) 

While Equation ( |11| ) makes the functional dependence of the solution more complicated than in 
Equation (12), it is useful for doing long time asymptotics, i.e. taking the limit as t 


oo. 


Equation (12) is useful for checking that the expression satisfies the equation. While the integrands 


of these expressions are different, the integrals are equal and thus one may switch between them 
whenever convenient. 

In Section [5j we show how it is possible to obtain expressions for all six spectral functions gj, 
hj in the relevant domains by solving a linear system. Indeed, for any r G {1, 3, 5}, if k G D R (the 
closure of D^ R ), then a certain number, say m, of the global relation equations 0 are valid for k. 
We must supplement these equations with 6 — m interface conditions to obtain a solvable system. 
Given the coefficients of gj, hj in Q, it is clear that the determinant of the linear system must be 
a polynomial in k. The criteria of Propositions [l]-[3] identify the cases in which this determinant is 
not identically 0, that is the system is full rank. For such a full rank system, it is always possible 
to choose R > 0 sufficiently large that D R contains no zeros of the determinant, which is essential 
in the proof of Proposition |4j We denote this linear system by 


AX = Y + e ik3T y, 


where 


X = ( g 0 (ik 3 , T),gi(ik 3 , T),g 2 (ik 3 , T), h 0 (ik 3 , T ), h x (ik 3 , T),h 2 (ik 3 , T)) 


(13) 


(14) 


explicitly (i.e., '*(■), j = 1,2 and non-homogenous terms from the interface conditions) and y 

which includes unknown expressions (i.e., j = 1,2). 


The right-hand side of (13) is expressed as the sum of Y. which includes expressions that are known 


5 Results 


In each of the following propositions, we assume that the interface conditions are not such that the 
problem reduces to a pair of BVP. It is a matter of trivial linear algebra to determine whether any 
particular problem has this property, and its well-posedness and solution are then known Eiza. 

In the case a\ > 0 and a 2 < 0, the phase velocity for x < 0 is positive and the phase velocity 
for x > 0 is negative. Thus, information from the initial conditions propagates toward the interface 
as in Figure [5] In this case we expect the minimal number of interface conditions to be necessary 
for a well-posed problem. 


Proposition 1. Assume a\ > 0 and a 2 < 0. Equation (13) is solvable for X if and only if two 
interface conditions are given. These conditions must be of the form 


(3 u q w (0,t) + (3r 2 q i3 \0,t) + /3i 3 <?i 3) (M) + /3i 4 <? (2) (0, f) +/3 15 q^\o,t) +Pwq^J{0,t) = hit), 


(15a) 

/3 2 ig (1) (o, t) + /3 2 2Q'i 1) (0, t) + faqW (0, t) + /W 2) (0, t) + /3 25 qM (0, t) + faq® (0, t) = f 2 (t). 

(15b) 
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t 



Figure 5: Information from the initial conditions q^\x, 0) and q ly ‘ 2 \x, 0) propagates toward the 
interface. 


The solution to (13) is full rank, that is, solvable for X, whenever at least one of the following 
holds 

1. j3 14/321 ~f~ /3ll/?24; 

2 . CTi(/3i5/?21 — /?11/^25) Y ^2(^12/324 — /3l4/322), 

3. erf (/3 i 6 /3 2 1 — P 11 P 2 Y) + eriCT2(/3l5/3 2 2 — ^ 12 / 325 ) + ^2 (/3l4/3 2 4 — /3l3^24) 7 ^ 0) 

-{• CTi(/3i 6 /3 2 2 — /3i2/3 2 Q) 7^ cr 2 (/3i3/3 2 5 — /3i5/323), 

h- ^16/323 Y /3l3/326- 

Proof of Proposition^ 1} In the case (Ti > 0 and 02 < 0, the second integrals of both ( |10a ) and (10b) 
can be deformed from • dA to — f dD (s) • dA. We rewrite the global relations ([7]) as 




o' 1 


f — ,T) - ql 1] ( — ) = a 3 g 2 (ik 3 ,T) + iaka(gi(ik 3 ,T) - (ak) z aig 0 (ik 3 ,T), (16a) 


aA 


o-i 


s* fc3t g( 2 ) ,T^ — q^ —'j = — a 3 h 2 (ik 3 ,T) — iaka^h^ik 3 ,T) + (ak) 2 a2ho(ik 3 ,T), (16b) 


ifc 3 t-(i) f ^ 


e q 


2 ; 


^,T ) - ( W— ) = afg 2 {ik 6 ,T) + ia 2 kafgi{ik 6 , T) - {a 2 k) z aigQ{ik 3 , T), (16c) 


a 2 A 


0 jfc 3 t~( 2 ) ( 


cr 2 


e**"*g w ( — ,T) -g® ( — ) = -^/ l2 (^ a ,T)-^aV^i(^ a ,T) + (a 2 A : )V2ho(^ a ,T), (16d) 


a 2 A’ 


0-2 


2 2 ; 


. 2 I\2. 


which are all valid for k £ z/ 5 \ Evaluating (15) for t = s, multiplying by e* fc3s , and integrating 
from 0 to t one obtains 

2 

hj(ik 3 ,T) + J2Pj+i,e+Mik 3 ,T) = f j+1 (ik 3 ,T), j £ {0,1,2}, 


e=o 


where 


.fj(oJ,T) = / e w 7i(s)ds, je {0,1,2}, 


10 










which is valid for k £ D^' 1 (the closure of D^). 

In order to solve the full 6x6 system it is clear we must impose two “interface conditions,” since 
the global relations (16) provide exactly four of the necessary six equations. If there is one boundary 
condition relating q^> and and their spatial derivatives then one can solve the problem on the 
left (right) and use the solution and remaining interface conditions to solve the problem on the 
right (left). The half-line problem is well posed [8j[22] and its solution will not be considered here. 
Hence, “interface conditions” of the type (15) are all we need to consider. 

/ c\ 

The above argument only fails if det(*A) = 0 since all singularities are outside D R . Examining 
det(„4.) = 0 one obtains a polynomial in k. Since we need this to hold for all k, we consider the 
coefficients of each power of k. Requiring at least one coefficient to be nonzero gives the conditions 
stated in Proposition ([!]). Q 


In the case ay > 0 and 02 > 0, the phase velocity for x < 0 and x > 0 is positive. Thus, 
information from q^\x) propagates toward the interface but information from q^\x) propagates 
away from the interface as in Figure [6j Hence, we expect that more interface conditions are 
necessary for a well-posed problem than in the case when oq > 0 and 02 < 0 as in Proposition [TJ 
Notice that the case of oq < 0 and 02 < 0 could be considered in this case by letting x —> —x. 
Hence, we consider only the case where cxi > 0 and 02 > 0. 


t 



Figure 6: Information from the initial condition q^\x, 0) propagates toward the interface while 
information from q^ 2 )(x,0) propagates away from the interface. 


only if three interface conditions of the following form are given. 


) is solvable for 

X if and 


(17a) 

/2W, 

(17b) 


(17c) 


The solution to ( |13[ ) is solvable for X whenever one or more of the following is satisfied: 

1. hi ± 0 , 

2 . oq/3 2 i + 0-2,032 / 0, 

3 . af /3 11 + 0-10-2/322 + 02#33 / 0, 
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4 - &iPi 2 + V2P23 / 0 , 

5. /3i 3 / 0. 


Remark. It may be possible to rewrite the interface conditions so that one is a boundary 
condition for c/ 2 - 1 and still have an interface problem. However, a single boundary condition for q W 
or a pair of boundary conditions for q^ 2 ' 1 implies that the problem separates into a pair of BVPs. 


Proof of Proposition In the case o\ > 0 and a 2 > 0, the second integrals of (10a) can be deformed 
from fff^-dk to — f ap ( 5 ) -d k. The second integral of (10b) can be deformed from ■ dk to 

R 

f 9D ( 1 ) • dk + f dD ( 3 ) • dk. We rewrite the global relations for each r E {1, 3, 5} as 


e ikH qW 




— ) = afg 2 (ik 6 ,T) + ika r cr(gi(ik 6 , T) - (ka r ) 2 a 1 g 0 (ik 6 ,T), 
Tl 


e ikH q W 


e ikh q( 2) 


a r+2 k 


v 1 


T ~q y 0 


(i) 


a 


' +2 k 


o-i 


a r+1 k 


02 


afg 2 (ik 3 ,T) + ika r+2 cr\g\(ik 3 , T) — (ka r+2 ) 2 aigo(ik 3 ,T), 
(2) / a r+1 JA _ 


T )~i°) 


which are all valid for k E D 
from 0 to t one obtains 


(18a) 


(18b) 


(18c) 


- a 3 2 h 2 (ik 3 ,T ) - ika r+1 a 2 hi(ik 3 , T) + (ka r+1 ) 2 a 2 h 0 {ik 3 ,T), 

Evaluating (|17[) for t = s, multiplying by e* fc3s , and integrating 


hj(ik 3 ,T) + /3 j+1/+ igc(ik 3 ,T) = f j+1 (ik 3 ,T), j E {0,1,2}, 


e=o 


where 


fj(u,T)= / e“ s f,(s) ds. j'E {0,1,2}, 
■Jo 


—(t*) 

which is valid for k E D 

In order to solve the full 6x6 system it is clear we must impose three “interface conditions”, 
since (18) provides just three equations. We must now examine the cases where one or more of 
these conditions decouples into a boundary condition on either q W or q( 2 \ If there is one boundary 
condition relating q l 1 ) and its spatial derivatives, then one can solve the problem on the left and 
use the solution and remaining interface conditions to solve the problem on the right. Solving the 
half-line problem is well posed [8, 22] and is not considered here. Hence, “interface conditions” of 
the type © are all we need to consider. 

Examining det(^l) = 0 in this case, one obtains a polynomial in k. Since we need this to hold for 
all k, we consider the coefficients of each power of k. Since we want conditions on det(yl) Owe need 
at least one of the coefficients to be nonzero. This gives the conditions stated in Proposition <§• 

□ 
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In the case o\ < 0 and cr 2 > 0, the phase velocity for x < 0 is negative and the phase velocity for 
x > 0 is positive. Thus, information from the initial conditions propagates away from the interface 
as in Figure [7j Hence, we expect that more interface conditions are necessary for a well-posed 
problem than in the previous cases. 

t 



Figure 7: Information from the initial conditions and 0) propagates away from the 

interface. 


Proposition 3. Assume o\ < 0 and cr 2 > 0. Equation (13) is full rank if and only if four interface 
conditions are given. These conditions must be of the form 


g (1) (0,t)+ $i4<? (2) (0, t) + $i 5 <t£ 2) (0, t) + $1 6 g£ } (0, t) 
4 1} (0, t) + $ 24 <? (2) (0, t) + $25<?i 2) (0, t) + $26?£z (0, t) 
QxJ (0, t) + $34<? (2) (0, t) + /?35^ 2) (0, t) + /?36gi 2) (0, t) 

$44 g (2) (0, t) + /?454 2) (0. *) + $464* (0, t) 



(19a) 

/ 2 W. 

(19b) 


(19c) 

f A (t)■ 

(19d) 


The solution to (13) is full rank whenever one or more of the following is satisfied: 
T P35 Pu / $34$45> 

2 . CT 1 ($34$46 — /336/?44) 7 ^ 02($24$45 “ $25 $44), 

3. 0 2 ($35$46 — /336/345) + 0102($26/344 — $24$46) + 0 2 ($14$45 _ $15$ 44 ) 7 ^ 0; 

0l($26$45 ~ $25$46) 7 ^ °2($16$44 — $14$46$ 

5. $16$45 / $15$46 • 


Remark. As four interface conditions are required, it must be possible to write (at least) 
two as boundary conditions. If there are two boundary conditions for either q^ or q^ 2 \ then the 
problem separates into a pair of BVP, so we only consider the case where there is precisely one 
boundary condition for each of q^ and However, for the purposes of stating the result, it is 
more convenient to write the conditions in the form (19). 
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Proof of Proposition^ In the case <t\ < 0 and a 2 > 0, the second integrals of both (10a) and (10b) 
can be deformed from ■ d k to J dD (i) ■ d k + f gD ( 3 ) • d k. The appropriate global relations can be 

rewritten for r G {1,3} as 


;T 

<J 1 J \ O'! 


e ik 3 t~(i) ( « r+lfc rr\ -(l) (ot r+1 k 


e lkH qW ( 


(20a) 


afg 2 (ik 3 ,T) + ika r+1 afgi(ik 3 , T) — (ka r+1 ) 2 crigo(ik 3 , T), 


a r+1 k \ 42) (a r+l k 


o-2 J (20b) 

- c r 3 2 h 2 {ik 3 , T) - ika r+1 a%hi{ik 3 , T) + (fca r+1 )V 2 /i 0 (iA; 3 , T), 

which are all valid for G D^\ Evaluating ( |19[ ) for t = s, multiplying by e lk3s , and integrating 
from 0 to t one obtains 

g 0 (ik 3 ,T) + p 15 h! (ik 3 , T) + p 16 h 2 {ik 3 , T) = h(ik 3 , T), 
gi(ik 3 ,T) + PMik 3 , T ) + p 26 h 2 (ik 3 , T) = f 2 (ik 3 , T), 
g 2 (ik 3 ,T) + fahi{ik 3 , T) + fo 6 h 2 (ik 3 , T) = f 3 (ik 3 , T ), 
h 0 (ik 3 , T) + PMik 3 , T) + P4sh 2 (ik 3 , T) = / 4 (ifc 3 , T), 


where 


-(?*) 

which is valid for k & D 


f j (u,T)= / e“ 
■Jo 


S fj(s)ds, je{ 1,2,3,4}, 


In order to solve the full 6x6 system (13) we must impose four “interface conditions,” since (20) 


gives only two equations. We need to examine the cases where one or more of these conditions 
decouples into a boundary condition on either qd') or . Using elementary linear algebra it is 
clear that at least one of these conditions must be a boundary condition. If there are two boundary 
conditions relating q^ ( q ( 2 )) and its spatial derivatives then one can solve the problem on the left 
(right) and use the solution and remaining interface conditions to solve the problem on the right 
(left). Solving the half-line problem is well posed (8)[22] and will not be considered here. 


There must be precisely one boundary condition of the form (19d). It is then elementary linear 


algebra to see that the remaining three conditions can be written as interface conditions in the 


form (19a)-(19c) 


Examining det(^4) = 0 one obtains a polynomial in k. We need this condition to hold for all k 
and we consider the coefficients of each power of k. Since we want conditions on det(^4) ^ 0 we need 
at least one of the coefficients to be nonzero. This gives the conditions stated in Proposition ([3]). □ 


Proposition 4. Assume A in (13) is full rank. A solution to 0 is given by ( [l2| ) where gj(ik 3 ,T) 
and hj(ik 3 ,T) for j = 0,1, 2 are the solution to the linear system AX = Y where A, X, and Y are 
given in (14) and the surrounding paragraph. 


Proof of Proposition [^} Consider Aj, which is the matrix A with the j th column replaced by e lk 1 y 


We solve AX = e lk3T y using Cramer’s Rule |5]. If we show that the contribution to the solution 
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from e lk3T y is zero, then we have proved the proposition. The terms we are concerned with 
from (12) are 


1/ 

2vr y r( i) 


x—ik 3 t 


(afg 2 (ik 3 ,T) + ikaigi(ik 3 , T) - k 2 g 0 (ik 3 ,T)) d k. 


and 


1 


i—x—ik 3 t 


e 


2vr J r ( 2 ) 

Using Cramer’s Rule and our factorization these become 


(cr|/r2(ifc 3 , T) + ika2hi(ik 3 , T) — k 2 ho(ik 3 ,T)) dfc. 


/ e 

J rt 1 ) 


1 

27T 
1 

2vr 7r( 2 ) 


i*+ifc3 (T _ t) 




2 det(^43) | det(yl2) 2 det(Ai)^ 


e CT 2 


^2 


det(y4) 

2 det(^.6) 

det(^4) 


+ ika\ 


det(„4) 


det(vA) y 


dA;. 




det(yl) 


det(A) y 


(21a) 


(21b) 


We would like to show these integrand terms are analytic and decay for large k inside the 
domains around which they are integrated. Note that det(A) ^ 0 since (13) is full rank. 

Case 1. (j i < 0, (72 > 0: For g\ < 0 and 02 > 0, = dD ^ U dD R \ Using the form of y 


in this case each term of the integrand in (21a) is of the form 

[ e ^+ ik 

J r(D 


1 

2n 

1 

2n 


j 

J rW 


i^+ik 3 (T-t) 
ik 3 (T-t)+ — 


m {1) 


a r+1 k 


,T ) +c 2 {k)q {2) 


a r+1 k 


02 


,T 


dk 


ci(k) 


7"+l 


q (1) (y,T)e ct i 


1 

+ 


2tt J p(i) 


«fc 3 (T—1) + — 
e v f "ic 2 


w (I 


,r+l. 


q (2 \y,T)e ^ 


dA; d y 
dk d y, 


where r 6 {1,3} depending on which region one is integrating around (D R ^ or D^) and ci(k) 
and C 2 (k ) involve the constants which are 0(1) as k —>• 00 from within r^ 1 ) and analytic 

for all k £ D R . For k —> 00 with kGD the expression inside the parenthesis decays by 
the Riemann-Lebesgue Lemma. Thus, by Jordan’s Lemma, these integrals along a closed, 
bounded curve in the complex k plane vanish for x < 0. In particular we consider the closed 
curves £ ■’) = C D (i) U£ C (i) and £^ = £n( 3 ) L)£ C ( 3 ) where C D u) = dD^ R n {k : \k\ < C} and 
£q(j) = {A: £ D^ R : \k\ = C}, see Figure 8 

Since the integrals along £ C (i) and £q(z) vanish for large C, the integrals must vanish since the 


contour C D {i) becomes dD R as C —> 00 . The same argument holds for £® and dD R . The 
uniform decay of the expressions in parentheses for large k is exactly the condition required 
for the integral to vanish using Jordan’s Lemma. Hence, (21a) is zero. A similar argument 
holds for (|21b|). 


Case 2. cji > 0,ct2 > 0: For o\ > 0 and (J 2 > 0, = dD R ^ and = dD^ R + dD R \ Using 
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Figure 8: The contours C D (i) and C D (z) are shown as red solid lines and the contours C C ( i) and 
C C ( 3 ) are shown as red dashed lines. The contour C D ( 5 ) is shown as a green solid line and the 
contour 5 ) is shown as a green dashed line. An application of Cauchy’s Integral Theorem using 
these contours allows elimination of the contribution of q^(-,t) and from the integral 

expressions (21). 


the form of y in this case each piece of the integrand in (21a) is of the form 
27r J r(i) 


i^+ik 3 {T-t) 


(k)qW 


crk 


1 


+ — / e" CT i 

Jr 


i^+ik 3 (T-t) 


2-k J r (i) 

ik 3 (T-t.)+i^ 


cr 1 
k 


1 T ) + C2 




c 3 m {2) ( —,t) dk 


+ 


-f ‘ 

27r J p(i) 

-/ 

2 -k J T ( 1 ) 

L 


Cl (k) 


ik 3 (T-t)+i^ 

e y ; CT i c 2 (fc) 


q^(y,T)e^j dkdy 

'‘O -ikocy) \ 

q (1 \y,T)e CT i I dfcdy 


+ ^ / e 

27T J r (i) 


ik 3 {T-t)+ 


ikx / / 

171 c 3 (fc) 


—00 

00 


— iky 


q {2) (y,T)e ^ ) dkdy, 


dk 


where c±(k), C 2 (k), and c 3 (k) are 0(1), analytic in D^ 5 \ and involve the constants /3j^. For 

k —> 00 with k G D the expression inside the parentheses decays by the Riemann-Lebesgue 
Lemma. Thus, we can apply Jordan’s Lemma and Cauchy’s Theorem as in Case 1, using 
curves shown in Figure [8j Hence, (21a) is zero. Again, a similar argument holds for (21b). 

Case 3 . o\ > 0,02 < 0 : For o\ > 0 and <72 < 0 , r^ 1 \r^ 2) = Using the form of y in this 


16 














case each piece of the integrand in (| 21 a[) is of the form 


1 

27T 


J e l 'T 1 +lk ( T d (fc)g (1 ) + c 2 (k)q (2) 






-/ 

Jr(i) 


ik 3 (T-t)+i 


ici(A:)( / q (1 \y,T)e CT i I d/cdy 


— ikay 


+ 2^ 


+ 2^ 


hL e ‘ 

IL/ 


f°° -ikay 


k ( T ^ + <o C 2 (fc) ( / q^iyij, T)e CT 2 ) dfedy 


fc (T * )+ 171 c 3 (fe) ( / q (1 \y,T)e CT i ) d&dy 


+ i /re, 6 


jfc 3 (T-t)+^ 


r —oo 
poo 


i C 4 (fc) [ / q ( ' 2 \y,T)e CT 2 ) dAdy 


—ika z y 


where c±(h), 02 (h), 03 (h) and 04 (A) are 0 ( 1 ), involve the constants /3j£, and are analytic for 

h E D^. For k 00 with A: E -D*' 3 '* the expressions inside the parentheses decays by the 
Riemann-Lebesgue Lemma. As before we apply Jordan’s Lemma to the appropriate curves in 
Figure[ 8 ]and use Cauchy’s Theorem. Thus, (21a) is zero. A similar argument holds for (21b). 


□ 


Remark. We show in Proposition [4] that ( | 11 [ ) is a solution of problem ([Tj) . It remains to show 
that this solution is unique in order to establish well posedness. In an attempt to show uniqueness, 
we assume there exist two solutions to 0- Let u(x, t) be their difference. Then u(x, t) satisfies 0 
with homogenous initial and interface conditions. A standard energy argument shows 


d t 


f 


d n 


dx ri 


d r 


u(x,t)\ dx = ai \2^—u (1) (0,t) 


dx r ‘ 


d n+2 .s / d n+1 


9 


- <T2 l 2 ^“ (2)(0 ’ t> 


dx n+2 

Qn+2 


U {1) (0,t) 


dx n+2 


u 


( 2 ) 


(0,0-( 


V dx n+1 

/ Qn +1 


u {1) (0,t) 


dx n+1 


u {2) (0,t) 


( 22 ) 


for any nonnegative integer n. If the interface conditions given are such that the right-hand side 
of ( 22 ) is always negative then, because u(a;, 0 ) = 0 and the left-hand side of ( 22 ) is always non¬ 
negative, we have u(x,t) = 0. Thus, one suitable choice of interface conditions is those that satisfy 
this relationship. For various signs of a\, 02 , which we consider here, it is not clear how to establish 
that the solution 0D> is unique in general. 


6 Examples 

In this section we give solutions to 0 for different signs of (J\ and 02 with “canonical interface 
conditions.” That is, we prescribe that the function and its first N spatial derivatives are continuous 
across the interface where 1 < N < 3 depends on the signs of a\ and oi- 
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Example 1. a\ < 0,a 2 > 0: This example requires four interface conditions. We impose that the 
function, as well as its first, second, and third derivatives are continuous across the boundary. 


qW(0,t)=qW(0,t), 

4 1} (° ,*) = 4 2) (m), 

= q£J(0,t), 

4**( ( m) = 4*Um), 

The first three conditions can be imposed directly. The condition on the third spatial deriva¬ 
tive can be imposed by applying the equation and integrating in t to give ([6]). Applying the 
t transform to (Jb]) we have 

\g 0 (ik\T) h 0 (ik 3 ,T ) = ClkT ~ l f44 1} ( 0) - 4^ 2) (°)) • ( 23 ) 

df (7% IK 6 \df (TJ J 

Using elementary row operations, we have, in the notation of Proposition [3j fi(t) = f 2 (t) = 
fo(t) = 0 , u(t) = e ’44 1 (^ 4 1 } ( 0 ) - ^ 4 2) (°)); #25 = /%6 = -1, and the remaining = 0. 
Using these interface conditions and solving (JT3]) , Equation (|TTj) becomes 


q {1 \x,t) = 


1 

27T 

+ . 
+ 

+ 


a 2 cri - <72 l) /a 2 fc 


aD (i) 27ra 2 cJi(cJi - <r 2 ) 


e CT i 


% 


cr 1 

ikx 


dk 


i{aX) -_ 1 j e 1^ g (i) (o) dk 


< dD W 2ira 2 ka 1 cr 2 (<7i - cr 2 ) 

IX 


L 


a;(a -!) .^-.^(2 ){*±\ ik 


SD ( b x) 27ra 2 cr|(cJi - cr 2 ) 


a 2 k 


cr 2 

ikx 


+ 


L 


an _ f Af d * 


+ 


+ 


9D ( d 3) 27rfeacJi(fJi - <72) 

i(a 2 - 1) ikH 


\ CT\ J 

(1) 


'dDf 2t rk{af - a f) 


(e"** * - l)eU g ^ (0) dfe 


d(« -1) ik 


‘dDW 27raa|((7i - a 2 ) 
f icr 2 {a 2 - 1) ifc s t 


\cr 2 J 


( e _ Jfc . t _ 1)e ^ g ( 2 ) (o) dfc 


IdD™ 2vrcr|aA:(c7i - cr 2 ) 
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« (2) 0M) = £ 


1 


ikx—LU2t ~(2) / 


2 7T 

+ 


e"™' "" 6 (/c) dA 


cr 2 (a 2 - 1 ) ihsi-ik 3 t„(i) ( « 2 A 


dD™ 2ira 2 crl(cr 1 - a 2 ) 


e CT 2 


96 


°i 


d A 


+ / *(«-l)(g2-qgi + gq 2 ) , ifc3t _ lk ^Ji) d ,. 

+ La. WfarJ^ - „ 2 ) (C 1)6 ®> (0) 


L 

L 


a 2 a 2 - ai i^-ik 3 t„( 2 ) / a 2 A’ 


9D a) 27ra 2 fj 2 ((Ti - cr 2 ) 


e a 2 


% 


02 


dA 


i(a — l)(acri — acr 2 — <r 2 ) , 


/9D 


C 1 ) 2na 2 ka2{o'i — a 2 ) 


(,e - l)eW q W (0) dA 


L 


cr|(a - 1) i£*~ik 3 t„( 1) (ak\ 


dD ( fi } 27racr 2 (cri - a 2 ) 


e CT 2 


% 


— dk 
O’! / 


Ua - 1)(lJ2 + a,7l) - l)e^(0) dk 


IdD™ 2'Kacr\k{ai - a 2 ) 


+ 


aa 2 - g i ik£- ik 3 t ^2) (ak\ 


l dD (3) 2naa 2 (a 1 - a 2 ) 


e CT 2 


% 


— dA 

cr 2 y 


L 


i(a-l)(^ + a» 1 ) (e _ a!St _ 1 )e t,f (0 ) di . 


az) (3) 27raAcr|(ai - cr 2 ) 


Remark. Combining ([ 6 ]) with the hrst interface condition of Example 1 and the mutual 
compatibility of c/q \x) and q^\x) with respect to the interface conditions, we find that 
q(0,t) is constant in time. That is, for all t > 0, = qj ) 1 ' 1 (0) = qjp (0) = q^\o,t). Note 

that this property of the interface problem is obtained without any reference to a solution 
method or any appeal to the first or second order interface conditions. The compatibility 
condition g^^O) = qj?\ 0) could also be used to simplify (23). We avoid this simplification 
because even without a smooth initial condition the formulae obtained in this paper are valid 
for t > 0 and the short time asymptotics can be analyzed as in [3J. 


Example 2. o\ > 0, c 2 > 0: This example requires three interface conditions. We impose that the 
function, as well as its first and second derivative are continuous across the boundary. That 
is, 


q^(0,t)=q^ 2 \0,t), 

4 1 } (0,0 = 4 2 ) (0,f), (24) 

4z(M = 4 2 J( 0 ,t). 


In the notation of Proposition [2] fi(t) = f 2 (t) = fs(t) = 0, fin = fi 22 = fi 33 = — 1, and 
the remaining fij^ = 0. Using the interface conditions (24) and solving (13), Equation © 


becomes 


19 













i roc 

q { 1 ) (x,t) = -n] k 


+ 


+ 


(cri - <t 2 )(<ti + acr i + aa 2 ) .(i) (ak\ A1 

e ai % I I d«, 

V^l / 


' 9Z) (5) 27racri(cJi - a 0 2 )( 0 i + cr 2 + acr 2 ) 


(J 2 — (Tl 


' 9z? ( 5) 27racri(cJi + 0 2 + a0 2 ) 


i^-ik 3 t .(i) I a 2 k 

a 1 9o I - 




dfc. 


.SCTj 5 


Jsd (5) 27T(T 2 (<7i - acr 2 )(cri + 0 2 + acr 2 ) 

R 

1 /*°° 

g (2) {x,t) = — / e ikx ~ U2t qQ {k) dfc 
2 tt 7_oo 

f er 2 {cr\ + cricr 2 - 0 2 ) 


i^-ik 3 t 4 2 ) 

e CT1 9 0 


IdD™ 2 ncrf{crf + a(l + a)cri(T 2 - 0 ^) 


i^~ik 3 t .(i) 
e CT 2 gg 


+ 


02(02(01 + (T 2 ) + Qf((TJ + 0-3)) !§M~ik 3 t .(1) 

e a 2 g 0 


+ 


/ai)W 27rafjf (of + a(l + a)critj 2 — 0 2 ) 
I" erf + (1 + a)0i0 2 - «0 2 

/ 9D U) 27ra0 2 (0 2 + a(l + a)aia 2 - erf) 
f ero(crf + er\cr 2 — erf) 


— ) d k, 
CT 2, 


—) dfc 

ff l/ 

d* 


'9D ( „ 3) 27raf(aerf(1 + a) + 02(01 + 02)) 


^~ik 3 t 42 ) 
e ^2 g^ 

— — *fc 3 t 41) 

e CT 2 g^j 


a 2 fe 

02 


dfc 


+ 


L 


0 2 (a 0 i( 0 2 - 01 ) + 0 2 (01 + 0 2 )) i±e-ik 3 t ,(l) 

—--------—e CT 2 q y 0 > 


dD ( 3 ) 2iracrf(acrf(l + a) + 0 2 ( 0 i + 0 2 )) 
f (1 + a) 0 2 + 0 i 0 2 + acr 2 

IdD^ 2naa 2 (aaf (1 + a) + 0 2 ( 0 i + 0 2 )) 


-) 4k 
<Tl) 

a 2 k 

01 


d k 


i^-ik 3 t 42 ) 
e CT 2 gg 


ik . 

02 / 


Example 3. 01 > 0, 0 2 < 0: This example requires two interface conditions. We impose that the 
function and its first derivative are continuous across the boundary. That is, 


q( 1 \ 0 ,t) = qW( 0 ,t), 

Qx\o,t) = qfX o,t). 


(25) 


In the notation of Proposition [I] fi(t) = / 2 (f) = 0, /J 15 = (3 2 5 = — 1, and the remaining 
f3j t e = 0. Using the interface conditions ( J25| ) and solving (13), Equation © becomes 
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pw> -1 £«^*4»wd* + i - 


/ 9d ( 5) 27raai((Ti + (J 2 ) 


aA;\ 


dfc 


+ 


+ 


02 + acr 2 - cti — —ifc 3 t ,(i) / a 2 A; 
--r a i q Q ’ I - 


'ao (5) 27racri(cri + 0 2 ) 




dfc 


cji(2 + a) — —ifc 3 t^( 2 ) /afc\ 


' 9z >( 5) 27ra<r 2 (cri + 0 2 ) 


e CT i 


96 


UJ* 


(Ji(2 + a) i^-ik 3 t 42 ) ( OL 2 k 


d D ( 5 ) 27racr 2 (cri + 0 2 ) 


e CT i 


96 


R 

J /*CJO 


5<2|(I ’ () = 2^/_ 


o ikx—0J2t 2) 


(/q ' (A?) dfc — 


— dfc, 

(j 2 


<t 2 (2 + a) ih^—ik 3 t ^(1) /aA;\ 
-e CT 2 % — 


'ao (5) 27ra(Ji((Ti + 0 2 ) 

R 


— dk 

) 


+ 


cr 2 (2 + a) — —jfc 3 tui) /a 2 A; 

- 

02 + acr 2 - 01 ih 2 L-ik 3 t 42 ) f ak\ 


'dD (5) 27ra0i(ai + a 2 ) 


dA; 


/ 9d ( 5) 27rao- 2 (0i + a 2 ) 


e CT 2 


96 


UJ* 


01 + «01 - 0 2 ih 2 L-ik 3 t 42 ) 
e CT 2 


/ 9d ( 5) 27ra0 2 (0! + 0 2 ) 

R 


a 2 k 

&2 


dk. 
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